Controllable manipulation and entanglement of macroscopic quantum states in 

coupled charge qubits 
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We present an experimentally implementable method to couple Josephson charge qubits and 
to generate and detect macroscopic entangled states. A large-junction superconducting quantum 
interference device is used in the qubit circuit for both coupling qubits and implementing the readout. 
Also, we explicitly show how to achieve a microwave-assisted macroscopic entanglement in the 
coupled-qubit system. 



I. INTRODUCTION 

Quantum-mechanical systems can exploit the funda- 
mental properties of superposition and entanglement to 
process information in an efficient and powerful way that 
no classical device can do. Recently, Josephson-junction 
circuits have received renewed attention because these 
may be used as qubits in a quantum computer. 1 Based 
on the charge and phase degrees of freedom in Josephson- 
junction devices, charge^ and phase qubits^^ have 
been developed. Also, a type of solid-state qubit can be 
realized in a large-area current-biased Josephson junc- 
tion^ 

Experimentally, coherent oscillations were demon- 
strated in a Josephson charge qubit prepared in a super- 
position of two charge states^ More recent experimental 
measurements^ showed that the charge qubit at suitable 
working points can have a sufficiently high quality of co- 
herence (Q v w 2.5 x 10 4 ), corresponding to a decoherence 
time T v w 500 ns. Current-biased Josephson junctions 
can also have long decoherence tiniest and Q v can reach 
10 4 . These exciting experimental advancements demon- 
strate the potential of Josephson qubits for manufac- 
turing macroscopic quantum- mechanical machines. To- 
wards the practical implementation of a solid-state quan- 
tum computer, the next important step would be the cou- 
pling of two qubits and then scaling up the architecture 
to many qubits. 

In this work, we present an experimentally imple- 
mentable method to couple two Josephson charge qubits 
and to generate and detect macroscopic quantum entan- 
gled states in this charge-qubit system. Motivated by 
very recent experimental results^ we employ a super- 
conducting quantum interference device (SQUID) with 
two large Josephson junctions to implement the readout. 
The generation of the macroscopic entanglement is as- 
sisted by applying a microwave field to each charge qubit. 
The key advantage of our design is that the SQUID can 
also produce an experimentally feasible and controllable 
coupling between the two charge qubits. As verified in 
a single qubit, - the coupled charge qubits may be well 
decoupled from the readout system when the measure- 
ment is not implemented. Moreover, our design can be 



readily extended to coupled multiple 10 qubits as well as 
any selected pairs (not necessarily neighbors). 

The paper is organized as follows. In Sec. II, the 
controllable coupling between two charge qubits is pro- 
posed using a large Josephson junction or a large-junction 
dc SQUID. Also, we demonstrate how this interbit cou- 
pling can be conveniently used to generate the controlled- 
phase-shift gate. In Sec. Ill, we study the microwave- 
assisted macroscopic quantum entanglement in the cou- 
pled charge qubits, where the microwave fields are cou- 
pled to the qubits via gate capacitances. Section IV 
focuses on the readout of the quantum states in the 
coupled-qubit system. Finally, the discussion and con- 
clusion are given in Sec. V. 



A. Other qubit coupling schemes 

A different type of interbit coupling from the one stud- 
ied here was proposed using the Coulomb interaction 
between charges on the islands of the charge qubitsiii 
As pointed out in Ref. the interbit coupling in this 
scheme is not switchable and also it is hard to make the 
system scalable because only neighboring qubits can be 
coupled. Implementations of quantum algorithms such 
as the Deutsch and Bernstein- Vazirani algorithms were 
studied using a system of Josephson charge qubits^ 
where it was proposed that the nearest-neighbor su- 
perconducting islands would be coupled by tunable dc 
SQUIDs. In Ref. 0, a pair of charge qubits were pro- 
posed to be capacitively coupled to a current-biased 
Josephson junction where, by varying the bias current, 
the junction can be tuned in and out of resonance with 
the qubits coupled to it. 

Another different type of interbit coupling was pro- 
posed^ in terms of the oscillator modes in an LC circuit. 
In contrast, we use a large junction or a large-j unction dc 
SQUID (but no LC circuit) to couple the charge qubits. 
In our scheme, both dc and ac supercurrents can flow 
through the charge-qubit circuit, while in Refs. and|3 
only ac supercurrents can flow through the circuit. These 
yield different interbit couplings (e.g., the a y a y typei^ as 
opposed to o x a x in our proposal). As revealed in Ref.Hol 
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and a capacitance C,/o. As in the single-qubit case^ close 
to the large Josephson junction, we also place a large 
capacitance Co in parallel with it, so that the effective 
charging energy of the large Josephson junction can be 
ignored (even though the capacitance of the large junc- 
tion might not be large enough). Moreover, we assume 
that the inductance of the qubit circuit (i.e., the two 
Cooper-pair boxes with the nearby junctions, and the 
superconducting lines connecting these two qubits with 
the large Josephson junction) is much smaller than the 
Josephson inductance of the large junction. The Hamil- 
tonian of the system can be written as 



FIG. 1: Schematic diagram of two charge qubits coupled by 
a large Josephson junction (denoted by a square with a "x" 
inside) of coupling energy Ejq and capacitance Cjo- To make 
the effective charging energy of the large Josephson junction 
as small as required, a large capacitance Co is placed close 
to and in parallel with it. Each filled circle denotes a su- 
perconducting island, the Cooper-pair box, which is biased 
by a voltage Vxi via the gate capacitance d and coupled to 
the bulk superconductors by two identical small Josephson 
junctions (each with a coupling energy Eji and a capacitance 
Co). Here the arrow near each Josephson junction denotes 
the chosen direction for the positive phase drop across the 
corresponding junction. 



the a x a x -type interbit coupling can be conveniently used 
to formulate an efficient quantum-computing scheme. 

Moreover, the calculated interbit-coupling terms in 
Refs. 0] and only apply to the case in which the fol- 
lowing two conditions are met: 

(i) The eigenfrequency ojlc of the LC circuit is much 
faster than the quantum manipulation frequencies. This 
condition limits the allowed number N of the qubits in 
the circuit because lolc scales with 1/^/N. In other 
words, this implies that the circuits in Refs.Q and0 are 
not really scalable. 

(ii) The phase conjugate to the total charge on the 
qubit capacitors fluctuates weakly. Our interbit-coupling 
approach discussed below is free from these two limita- 
tions. 



II. CONTROLLABLE COUPLING OF CHARGE 
QUBITS 

A. Coupling qubits with a large junction 

We first use a large Josephson junction to couple two 
charge qubits (see Fig. ^) . Each qubit is realized by a 
Cooper-pair box, where a superconducting island with 
excess charge Qi = 2ehi (i — 1, 2) is weakly coupled to 
the bulk superconductors via two identical small junc- 
tions (with Josephson coupling energy Eji and capaci- 
tance Cji) and biased by an applied voltage Vxi through 
a gate capacitance Cj. The large Josephson junction on 
the left has a coupling energy Ejq (much larger than Eji) 



11 = ^ [ E a(ni - n X iY - Eji(cos(j)iA + cos^ iB ) 

i=l 

~Ej cosj, (1) 



where 



E r 



2c 2 



Ci + 2Cji 



(2) 



is the charging energy of the superconducting island and 
nxi = CiVxi/2e is the reduced offset charge (in units 
of 2e) induced by the gate voltage. Flux quantization 
around loops containing the phase drops of the involved 
junctions gives the constraint: 



t>iA - <i>iB - 7 + 



2tt$ p 



0, » = 1,2, (3) 



which gives 
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(4) 



where the average phase drop 0, = \ {<t>iA+<t>iB) is canon- 
ically conjugate to the number, hi, of the excess Cooper 
pairs on the ith superconducting island: 



i = i,2. 



Here <piA and 4>iB (i = 1,2) are the phase drops across 
the small Josephson junctions above (A) and below (B) 
the ith Cooper-pair box. 

The Hamiltonian can be rewritten as 



H 



E 



E cl {hi - n X i) 2 - 2Eji cos 
- Ejq cos 7. 



$0 



X COS ( 



(5) 



The externally applied flux <i> e threads the area between 
the large Josephson junction and the left Cooper-pair 
box. It induces circulating supercurrents in the qubit 
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circuit. The total circulating supercurrent / has contri- 
butions from the two charge qubits: 



I = h+h 



(6) 



where 



'7T<£> e 1 \ 

I, = 21, , sin ( — -7) cos0j, (7) 

with I c i = TrEji/<S>Q. This total supercurrent flows 
through the large Josephson junction and it can also be 
written as 



/ = losing, 



(8) 



with I = 2irEj /$ . From Eqs. (6)-(8) it follows that 



Iq sin 7 = 2 sin 



$ 



- o7 ) (-^cl COS 01 + I c2 COS 2 ). 



(9) 

When the coupling energy Eji = ^ n I ci /TT of each Joseph- 
son junction connected to the charge box is much smaller 
than that of the large Josephson junction in the circuit, 
the phase drop 7 across the large junction will be small. 
Expanding the operator functions of 7 in Eq. JUJl into a 
series and retaining the terms up to second order of the 
parameters 



y*(<l), i = l,2, 



(10) 



we have 



: = - } sin ( -q-^-j (vi cos <t>i + cos 2 



$ ) \ 



T]i COS 01 + 1] 2 COS 0; 



(11) 



It is clear that the phase drop 7 across the large Joseph- 
son junction is controllable via the applied flux <f> e . 

For Hamiltonian (JSJ , we also expand the operator func- 
tions of 7 into a series and retain the terms up to second 
order of rji. Moreover, we consider the charging regime 
with E c i much larger than Eji. Also, we assume that the 
temperature is low enough (fc^T <C E ci ) and the super- 
conducting gap is larger than E c i, so that quasiparticlc 
tunneling is strongly suppressed. In this case, only the 
lowest two charge states are important for each qubit op- 
erating around the degeneracy point Vxi = (2rii+l)e/Ci. 
In the spin-i representation based on the charge states 
\rii) = I and |rij + 1) = | J,)j of each Cooper-pair box, 
the Hamiltonian of the system can be reduced to 



H 



with 



E [ e i<y*i) 



of-Ej^ 



x4^\ (12) 



£i{V X i) = -E c 



CiV; 



ivxi 



(2n, + 1) 



(13) 



and 



where 



E.i = E, <-<^ [ ) Li- 



& = l-^+3^)sin*( ^ 
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(14) 



(15) 



and i,j = 1,2 (i ^= j). The interbit coupling \ is given 
by 



X = Ljl c il c2 sin 



7T<I> C 



(16) 



where the large Josephson junction acts as an effective 
inductance of value 



L., 



$0 
2vr/ n 



(17) 



It is clear that the interbit coupling is switched off at 
$ e = 0. It is well known that a large Josephson junction 
can act as an inductance (e.g., Ref.[l|). Here we explicitly 
show a specific way that it can be used to couple qubits. 

Retaining up to second-order terms in the expansion 
parameters rji , the total circulating current / can be writ- 
ten as 

/ = 2sin I -j— J (I c i cos 0i + I c2 cos0 2 ) 



1 /2tt$ £ 
— sin — — 



(Id COS 01 + I c2 cos0 2 ) 2 - (18) 



In the spin-i representation, it is given by 
7 = sin(^) (Id^+I^) 



1 . /2tt$ £ 

sin 

47 V *o 



(l)+2) 



(19) 



which depends on the states of the charge-qubit system. 



B. Coupling qubits with a SQUID 

There are somewhat conflicting requirements imposed 
on this circuit. To obtain a large value for the effective 
Josephson inductance Lj = $o/27r/o, a relatively small 
Io is needed, so that a large interbit coupling can be 
achieved. However, when the large Josephson junction is 
also employed for a readout, it is desirable to use a large 
Iq. This permits a larger range of so that a higher 
resolution in distinguishing qubit states can be achieved 
in the quantum measurement based on the switching of 
the supercurrent through the large junction. 

These two opposite requirements can be conveniently 
solved if the leftmost large Josephson junction in Fig. 2] 
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FIG. 2: Schematic diagram of the coupled-qubit circuit with 
a biased-current source of impedance Z{u). The dc SQUID, 
with two junctions of large _E} , plays the role of both coupling 
the charge qubits and implementing the readout. Here the 
large capacitance Co placed close to and in parallel with the 
dc SQUID is included in the impedance Z(uj). 



is replaced by a symmetric dc SQUID with two suffi- 
ciently large junctions (see Fig. [2J. Instead of $ e inside 
the circuit loop between Ejq and the first qubit (as in 
Fig. 1), we now apply a flux $ s inside the large-junction 
dc SQUID loop (see Fig. 2). This SQUID can be used 
both for coupling the two charge qubits and implement- 
ing the readout. When the readout is not active (lb = 0), 
we can choose a suitable flux <J> S inside the SQUID loop 
to generate a larger interbit coupling. For It — 0, the 
reduced Hamiltonian of the coupled-qubit system and 
the total circulating current / have the same forms as 
in Eqs. ((T^|l and i|T$|) . but with $ e and I replaced by 
i$ s and 



In = 2/n cos 



(20) 



where Iq = 2irEj Q /& 



When the readout is active (see 
Sec. IV), <I> S is chosen as zero to obtain a larger effective 
Josephson coupling energy. 



C. Controlled-phase-shift gate 

When the system works at the degeneracy points with 
£i(Vxi) = 0, the Hamiltonian becomes 



H=-E J1 a^-E J2 a^- X a^a^. (21) 



For instance, when Eji > 0, i = 1, 2, its four eigenvalues 
are 



Eji + Ej2 - x, 
Eji - E,] 2 + x, 
Ej2 — Eji + x, 
—Eji — E J2 - x- 



The corresponding eigenstates are |ei,e2), \ei,g 2 ), 
\gi,e 2 ), and \gi,g 2 ), where 

|e i ) = ^(|T) i -|4) i ), 

lfli> = ^(lt> i + li>i). (23) 

Because they are also the eigenstates of the two uncou- 
pled charge qubits, when prepared initially at an eigen- 
state, the system does not evolve to an entangled state 
even in the presence of interbit coupling. As shown be- 
low, one can take advantage of this property to imple- 
ment the measurement. In addition, this property can 
be used to construct efficient conditional gates. For in- 
stance, if 

En = E.,2 = x, (24) 
the controlled-phase-shift (CPS) gate is given by 

C/cps(t) = e^ n U, (25) 

with 



U : 



-iHr/h 



exp 



{l ( X r/h)[a^+a^+a^a^}}, (26) 



at r = ttTi/Ax- This gate transforms the basis states 
|ei,e 2 ), \ei,g 2 ), \gi,e 2 ), and \gi,g 2 ) as 




/ 1 
10 
10 

\ -1 




(27) 



The generation of this conditional two-bit gate is effi- 
cient because the condition (24) can be realized in one 
step via changing the gate voltages Vx%, i = 1,2, and the 
flux <& s simultaneously. Also, the architecture is scalable 
because multiple charge qubits can be coupled by con- 
necting them in parallel with the large-junction SQUID. 
If the two Josephson junctions in each Cooper-pair box 
are replaced by small-junction dc SQUIDs, any selected 
pairs of charge qubits (not necessarily neighbors) can be 
coupled^ 



III. MICROWAVE-ASSISTED MACROSCOPIC 
ENTANGLEMENT 

When a microwave field is applied to the Josephson 
charge qubit, Rabi oscillations occur in the system^ 
These oscillations can also be demonstrated by coupling 
a quantum resonator to the charge qubit 1-^ Here we ap- 
ply the microwave field to the Cooper-pair box via the 
gate capacitance, as in Refs. andlliL but each charge 
qubit is driven by a different microwave fields In this 
situation, nxi in Eq. Q is replaced by 



(22) 



nxi + fiACi = n X i + 



2e 



(28) 



5 



Here di is the thickness of the gate capacitor and 

£aci = £\ia>i + £Ii<4 ( 29 ) 

is the microwave electric field in the gate capacitor of 
the ith Cooper-pair box, where cii is the annihilation op- 
erator of the microwave mode. Because the microwave 
wavelength is much larger than di , £ \i can be considered 
constant in the gate capacitor. In the charging regime, 
the Hamiltonian of the system (including the microwave 
fields) can be written as 



H 



Xi a. 



hu>\i aid 



i=l 



+af (Km + K*4)] - xa^af, (30) 



where 



Ki = 



EdCidi 
2e~ 



£\i 



(31) 



Here, we also consider the system working at the degen- 
eracy points Si(Vxi) = 0, i = 1,2. When hui\i « 2|J5ji| 
and under the rotating-wave approximation, the Hamil- 
tonian is cast to 



H = 



n- 



Eji °i + ^ w Ai aid} + (Ki\ei)(gi\ai + H.c.) 



(32) 



Without interbit coupling, each Josephson charge qubit 
exhibits Rabi oscillations between states |ej, Z<) and 
\gi, h+1), where \h) is a photon state with h photons. For 
the resonant case with hu>\i = 2\Eji\, the eigenvalues of 
each charge-qubit system are given by 



Si) 



where 



and 



E 0i = hoJ\i(k + 1), 



«i = T\K t \Vh + l 



(33) 



(34) 



(35) 



is the Rabi frequency. Though entanglement occurs be- 
tween each charge qubit and the nonclassical microwave 
field, the two qubits do not entangle with each other since 
the system evolves as 

\*®) = \Mt))\Mt))> (36) 

where 

\ipi(t)) = sin(f^)|e 4 , k) + cos(n z t)\gi, h + 1) (37) 

if the system is initially prepared at state \gi,g 2 , h+1, h+ 
1). However, in the presence of microwave fields, when 



the interbit coupling is switched on, the coupled-qubit 
system exhibits complicated quantum oscillations and it 
will evolve to the entangled state. For instance, in the 
resonant situation, the eigenvalues are given by 



£1,4 = Eoi + E 2 ± ftAi, 

£2,3 = Eqi + Eq2 ± hA2, 



where 



Ai, 2 = [(fix ± n 2 ) 2 + (x/Kf] 1/2 . 
The state of the coupled-qubit system evolves as 



(38) 



(39) 



|*(t)) = C 1 (t)\e 1 ,e 2 ,li,l 2 ) +C 2 (t)\e 1 ,g 2 ,h,l2 + 1) 
+C 3 (t)\g 1 ,e 2 ,h + l,l 2 ) 
+C 4 (t)\ 9l ,g 2 ,h + l,l 2 + l). (40) 

For the system prepared initially at \gi,g 2 , h + 1, Z2 + I}, 



Ci(t) = 


\{R2(t) 


-Ri{t) + 


-ib/h)[s 2 (t)-Si{t)]}, 


c 2 (t) = 






+ (fii-n a )5 2 (t)], 


c 3 (t) = 






- (ni-n 2 )52t)], 


C 4 (t) = 




+ R 2 (t) -+ 


- t (x/h)[s 1 (t) + s 2 (t)}}, 



where 



Ri(t) = cos(A i t), Si(t) = 



sin(A^) 
A,- ' 



(41) 



(42) 



For a two-level system interacting with a single-mode 
field, the Rabi oscillations can be explained using ei- 
ther quantum or semiclassical theory, where the single- 
mode field is described quantum mechanically or treated 
as a classical fields Here the quantum oscillations of 
coupled charge qubits (namely, the Rabi oscillations 
in coupled two-level systems) are studied using quan- 
tum theory, where the microwave field coupled to each 
qubit is quantized. This also applies to the classical- 
field case, in which the quantum oscillations are still de- 
scribed by Eq. gUJ), but |ei, e 2 , l\, l 2 ), \e\,g 2 ,hM + 1), 
\gi,e 2 ,h + l,h), and \g\,g 2 ,h + l,h + 1) are replaced 
by |ei,e 2 ), \e x ,g 2 ), \gi,e 2 ), and \gi,g 2 ). 

Figure shows the occupation probability |Ci(i)| 2 as 
a function of time t. For instance, when |Ci(t)| 2 ~ 1, 
both charge qubits are in their excited states. It can 
be seen that |Ci(i)| 2 looks very different when the in- 
terbit coupling is switched on or off. The macroscopic 
entanglement between the two coupled qubits can be ex- 
plicitly shown at f2i = Q 2 (= In this case, when 
tent = nirh/Wx, with n — 1,2,3,..., and 



w = [(2hn/ x ) 2 + i} 1/2 , 



(43) 
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FIG. 3: Occupation probability |Ci(t)| 2 as a function of time, 
(a) ft 2 = fii.x/ft = fli! (b) n 2 = 1.2«i, x A = Oi; (c) n 2 = 
Oi,x/^ = \/3«i/2; (d) n 2 = 1.2fii,x/7i = v^fii/2; (e) 2 = 
f2i,X = 0; (f) n 2 = 1.201, x = 0. The time is in units of Of 1 . 



Now, a term — <I>o/fc<5/27r, with 

E 



(47) 



should be added to the Hamiltonian QJ, where j is the 
average phase drop of the total qubit circuit and it can 
be written as 



7 2*o ' 



(48) 



with 7 = + 7h)- Here we set the flux <& s equal to 
zero to have a larger effective Josephson coupling energy. 
In the spin-i representation based on charge states, the 
Hamiltonian of the system is also reduced to Eq. JTJJ. 
The interbit coupling is here induced by the bias current 
and given by 



X = LjI cl I c2 sin 2 (7o/2), 
where the effective inductance is 



L, 



2ttIq cos 7q ' 



and 



7o = sin 1 (I b /I ), 



(49) 



(50) 



(51) 



|^(t)) becomes 

l*(*ent)) 



Ci(£ C nt)|ei, e 2 ,h, h) 

+ Ci (tent) | gi ,92,h + l,h+ 1) 



(44) 



where 



Cl(tent) 



[— cos(n7r) + exp(wwr/W)], 
[cos(n7r) + exp(in7r/TT r )]. 



(45) 



The peaks away from either zero or 1 shown in Fig. |3ja) 
correspond to this kind of entangled state. Furthermore, 
if suitable values of W are taken, the maximally en- 
tangled state with | C x | 2 = \Ci\ 2 = | can be derived. 
This state is a macroscopic Schrodinger-cat state of the 
two charge qubits. For instance, if hfl/x = V3/2, the 
coupled-qubit system evolves to the maximally entangled 
state at the times given by 



4X X) = (21 + l)nH/2x, 1 = 0,1,2,. 



(46) 



This entangled state corresponds to the half-probability 
peaks in Fig. Etc). 



IV. QUANTUM MEASUREMENT 

To implement a readout, we bias a current pulse to 
the qubit circuit (see Fig. EJ) , as in the single-qubit caseM> 



with Iq = 4irEj /&o, and h < Iq. The intrabit couplings 
are 



where 



with 



Ejt = Eji cos(7o/2)£i 



e i = l-a(r 7 ?+3r ? ?)sin 2 (7o/2), 



2 + cos 70 
8cos 3 7o 



(52) 



(53) 



(54) 



and i,j — 1,2 (i ^= j). The supercurrent through the 
SQUID, 

J sin 7 = lb - sin(7o/2)[/ c icr^ ) + I&cr^] 



4/ 



tan 7q [7^1 



l c2 + ZI c il c2 a x o x J, 



(55) 



has contributions from both the bias current and the cur- 
rent from the Josephson charge qubits. 

At the working points with Ei(Vxi) — 0, the eigen- 
states of the system are also \e\,e 2 ), \e\,g 2 ), \gi,e%), 
and |3i,32)- In Fig. we show the dependence of the 
supercurrents through the SQUID on the eigenstates of 
the charge-qubit system. The supercurrents through the 
SQUID increase with the bias current and the difference 
between the supercurrents at different (nondegenerate) 
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FIG. 4: Eigenstate dependence of the supercurrent through 
the SQUID as a function of the bias current h- Here, Eji — 
E J2 = \E%, jl> = | ei ,e 2 ), |2) = |e ljff2 ), |3> = \gi,e 2 ), and 
|4) = |ffi,£te). 

eigenstates widens. For the measurement setup shown 
in Fig. [3 the supercurrent through the SQUID is the 
largest at the eigenstate |ei,e2) and it first reaches the 
maximal value Iq (namely, the critical current) when 
the bias current lb approaches a value Isw near Iq. 
Around this value, the supercurrent through the SQUID 
switches, with a very large probability Pi, from the zero- 
voltage state to the dissipative nonzero-voltage state in 
the quasiparticle-current branch and the measurement 
on the voltage is carried out. However, due to envi- 
ronmental noise as well as thermal and quantum fluc- 
tuations, the switching actually occurs before the super- 
current through the SQUID reaches Iq. At lb ~ Isw, 
the supercurrents through the SQUID will also switch to 
the nonzero voltage state at other eigenstates, but the 
switching probabilities are small. In the ideal case, if the 
difference between the large switching probability Pi and 
the small ones is close to 1, then, in principle, a single- 
shot readout would be achieveable. As shown in Ref. 0, 
the Josephson-junction switching experiment can provide 
sufficient accuracy to discriminate the state |ei, e<i) from 
others. 

The operation and readout of the macroscopic en- 
tanglement of the coupled-qubit system can be imple- 
mented by simultaneously applying a pulsed microwave 
field (with the same duration r) to each charge qubit. 
The sequence would be: 

(i) before the microwave fields are applied, the flux $ s 
through the SQUID is set equal to zero and no interbit 
coupling exists; 

(ii) the flux $ s is switched on to a certain nonzero value 
exactly at the start of the microwave pulse and off at the 
end of the microwave pulse. Within the microwave pulse 
duration r, the evolution of the system is described by 
Eq. (U; 



(iii) a pulsed bias current lb is applied to perform a 
measurement after the microwave pulse. 

During the measurement, the quantum state of the 
charge-qubit system collapses to the eigenstate \ei,e2) 
with probability |Ci(r)| 2 . This probability is propor- 
tional to the switching probability Pi of the SQUID. Be- 
cause of relaxation, the envelope of the measured switch- 
ing probability Pi decays exponentially with time. This 
is used to obtain the relaxation timem& Ramsey fringes 
of the probability Pi can be used^ to determine the de- 
coherence time of the coupled-qubit system. For each 
given microwave pulse duration r, through repeated mea- 
surements, one can determine the occupation probabil- 
ity |Ci(t)| 2 and thus deduce the information about the 
macroscopic entanglement between the coupled charge 
qubits [see Figs. 3(a) and 3(c)]. 

V. DISCUSSION AND CONCLUSION 

Finally, we estimate some important parameters using 
available quantities for the single charge qubit. Here we 
consider the maximally entangled case shown in Fig.[3Jc), 
in which Oi = f?2 = and 



Taking 2tt/H, ~ 0.22 /its, as derived from the Rabi oscil- 
lation of the measured switchingprobabilityfS we have 
x/fi ~ 0.25 GHz. Reference |!J also gives Ej/h ~ 
16.5 GHz. Choosing 

E a J0 w hE Jl w 5P/, 

and using the relation for x, we obtain $ s 0.35$o. For 
f&s = 0, the expansion parameters are 

m = !°L « 0.05 
Jo 

for E S J0 w 5Eji. When $ s w 0.35$ , they become rfr w 
0.14. The results are sufficiently accurate when Eji and 
X are retained up to second- and higher-order terms in the 
expansion parameters rji. When Q s approaches $o/2, the 
interbit coupling strengthens. The reduced Hamiltonian 
of the system also has the same form as Eq. (|12|) . but 
higher-order terms in the expansion parameters should 
be included to obtain accurate results. 

Here we consider the charging regime with E C i ^> Eji 
in order to obtain analytical results. We expect that the 
interbit coupling can still be realized in the regime with 
E C i ~ Eji, i.e., the regime used by the Saclay group in 
the experiment on a single Josephson qubit ^ In this lat- 
ter regime, the results can only be obtained numerically, 
but a relatively long decoherence time would be expected 
for the coupled-qubit system to work at the degeneracy 
points because at these points the states are more stable 
against the variations of both the offset charges and the 
flux <I> e or $ s . 
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Very recently, quantum oscillations were experimen- 
tally observed in two coupled charge qubitsi^ Also, a 
novel method for the controllable coupling of charge 
qubits was proposed using a variable electrostatic trans- 
former. 19 In contrast with our interbit coupling scheme, 
these studies involve capacitively-coupled (as opposed to 
inductively-coupled) charge qubits. The main advantage 
of this inductive coupling among qubits is that it allows 
a controllable link between any selected qubits, not nec- 
essarily nearest neighbors. 

In conclusion, we employ a large-junction SQUID to 
couple Josephson charge qubits and implement a readout. 
This architecture is readily scalable to multiple qubits. 
When the system works at the degeneracy points, where 



the dephasing effects are suppressed, it is shown that the 
macroscopic entanglement can be generated with the as- 
sistance of microwave fields. Also, we show the quantum 
measurement of the macroscopic entanglement. 
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